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We study the one dimensional branching Brownian motion starting at the origin and investigate 
the correlation between the rightmost (Xmax > 0) and leftmost (Xmin < 0) visited sites up to time 
t. At each time step the existing particles in the system either diffuse (with diffusion constant D), 
die (with rate a) or split into two particles (with rate b). We focus on the regime b < a where 
these two extreme values Xmax and Xmin are strongly correlated. We show that at large time t, 
the joint probability distribution function (PDF) of the two extreme points becomes stationary 
P{X, Y,t ^ oo) —^ p(X, Y). Our exact results for p{X, Y) demonstrate that the correlation between 
Xmax and Xmin is nonzero, even in the stationary state. From this joint P DF, w e compute exactly the 
stationary PDF p((() of the (dimensionless) span (( = (Xmax — Xmin)/\/ D/b, which is the distance 
between the rightmost and leftmost visited sites. This span distribution is characterized by a linear 
behavior p(0 ~ i (1 + A) <(■ for small spans, with A = (^ — l). In the critical case (A = 0) 
this distribution has a non-trivial power law tail p(^) ~ 87r%/3/C^ for large spans. On the other 
hand, in the subcritical case (A > 0), we show that the span distribution decays exponentially 

as p(C) ~ (A^/2)('exp ((j for large spans, where A is a non-trivial function of A which we 

compute exactly. We show that these asymptotic behaviors carry the signatures of the correlation 
between Xmax and Xmin. Finally we verify our results via direct Monte Carlo simulations. 
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I. INTRODUCTION 


Branching Brownian motion (BBM) is a well-known 
model that finds applications in several areas of sci¬ 
ence including physics, mathematics and biology. BBM 
arises naturally in the context of systems where new par¬ 
ticles are generated at each time step such as models 
of evolution, epidemiology, population growth and nu¬ 
clear reactions, and now has a long history [ll-[l^. In 
addition, BBM has also been widely used in theoreti¬ 
cal physics where it has been studied in the context of 
reaction-diffusion models, disordered systems amongst 
others [H El- BBM is also an important model in 
probability theory as it combines the long-studied dif¬ 
fusive motion with the random branching mechanism of 
Galton-Watson trees da. In this paper we are interested 
in one dimensional BBM. The process begins with a sin¬ 
gle particle at the position a: = 0 at time t = 0. The 
dynamics proceeds in continuous time, where in a small 
time interval At, each particle splits into two independent 
particles with probability b At, dies with with probability 
a At, and with the remaining probability (1 — (a -I- 6) At) 
performs a Brownian motion on a line with a diffusion 
constant D. A realization of the dynamics of such a pro¬ 
cess is shown in Fig. [T] 

In a given realization of this BBM process, there are 
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FIG. 1: A realization of the dynamics of branching Brownian 
motion (left) in the supercritical case (right) in the critical 
case. The span of the process is defined as s = Xmax — Xmin, 
where Xmax and Xmin are the maximum and minimum dis¬ 
placements of the process up to a certain time t respectively. 


in general N{t) > 0 particles present in the system at a 
particular time t. The parameters b and a in this BBM 
model define three regimes with different properties. The 
number of particles N(t) is a random variable whose 
statistics depends on a and b. When the rate of birth 
is greater than the death rate (a < b), the supercritical 
phase, the process is explosive and the average number 
of particles in the system grows exponentially with time 
(X(t)) = exp((6 — a)t). In contrast, when the birth rate 
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is smaller than the death rate {b < a), the subcritical 
phase, the process eventually dies and, on an average, 
there are no particles present in the system as t — > cxd. 
At the critical point a = b, the system is characterized 
by a fluctuating number of particles with {N{t)) = 1 at 
all times t. 

If one takes a snapshot of the system at a given time t, 
the spatial positions of the existing particles happen to be 
strongly correlated, since the particles are linked by their 
common genealogy. One important object that has been 
extensively studied is the order statistics of these parti¬ 
cles, i.e., the statistics of the position Xk{t) of say the A:-th 
rightmost particle at time t, where the particle positions 
on the line are ordered as xi{t) > X 2 {t) > .... > a;Ar(t)(t) 
i ii, E E. Another related interesting quantity 
is the gap gk{t) = Xk(t) — Xk+i{t) between the k-th 
and {k + l)-th particle at time t. Most of these stud¬ 
ies have thus focused on extreme value questions at 
a given time t. However, there are other interesting 
extreme value observables that concern the history of 
the process over the entire time interval [0,t]. For in¬ 
stance, one can consider the global maximum, ATmax = 
maxo<.r<t [{a;i(r),a: 2 (T),X 3 (T)...,a::Ar(,-)(T)}] which rep¬ 
resents the maximum of all the particle positions up to 
time t. This has the simple interpretation as the max¬ 
imum displacement of the entire process up to time t 
(see Fig. [1]). This global maximum has appeared in 
a variety of applications including the spread of gene 
populations and the propagation of animal epidemics 
in two dimensions 0- Similarly the global minimum 
Amin = mino<r<t [{a:i(T),a:2(T),a:3(T)...,a;jv(r)(T)}] is 
another interesting quantity that, by symmetry, has the 
same marginal probability distribution function (PDF) 
^S Amax- 

The marginal PDF of Amax has been studied exten¬ 
sively for the supercritical , critical and the sub- 
critical phases [J, 0. While the marginal distributions 
of Amax, and hence that of —Xmin, are well studied, 
much less is known about the correlation between these 
two random variables. In this paper, we study the joint 
PDF of Amax and Amin- In the supercritical phase, this 
joint PDF is always time dependent, and is hard to com¬ 
pute analytically. However, in this case, Amax and Amin 
get separated from each other ballistically in time and 
hence become uncorrelated at late times. In contrast, 
in the critical and subcritical phases (b < a), we show 
that the joint PDF reaches a limiting stationary form 
at late times, which we compute analytically. More¬ 
over, for b < a, our exact results for the stationary joint 
PDF demonstrate that this correlation between Amax 
and Amin remains finite even in the stationary state. 

The joint PDF of Amax and Amin has the following in¬ 
teresting physical application. For instance, in the con¬ 
text of epidemic spreads, it is important to characterize 
the spatial extent over which the epidemic has propa¬ 
gated up to time t. This is clearly measured by the span 
s = Amax — Amin of the process up to time t (see Fig. [1]) 
[2l| - |^ . Evidently, to compute the distribution of s, we 



FIG. 2: Schematic representation of a trajectory of the BBM 
confined in the box [—Y, X]. Note that Amax and Amin denote 
respectively the maximum and the minimum of the process 
up to time t. The process starts with a single particle at the 
origin at time t = Q and hence Amax > 0 while Amin < 0. 

need to know the joint PDF of Amax and Amin- In this 
paper we also compute analytically the stationary PDF 
of the span s in the critical {b = a) and the subcritical 
(6 < a) cases. Our exact results demonstrate that the 
correlation between Amax and Amin is also manifest in 
the stationary span PDF. 

The rest of the paper is organized as follows. In sec¬ 
tion H, we define the model precisely and summarize our 
main results. In Section HI we derive an exact evolution 
equation for the joint distribution of Amax and Amin- In 
Section IV we derive the stationary joint PDF of Amax 
and Amin for the critical (6 = a) and the subcritical 
(6 < a) cases. In Section V we compute the stationary 
PDF of the span and extract its asymptotic behaviors 
analytically. In Section VI we compare our analytical 
predictions with Monte Carlo simulations. Finally, we 
conclude with a discussion in Section VH. Some details 
of computations are relegated to the appendices. 

II. THE MODEL AND A SUMMARY OF THE 
RESULTS 

The model and the observables. We consider the 
BBM on a line starting with a single particle at the ori¬ 
gin at time t = 0. The process evolves via the following 
continuous time dynamics. In a small time interval At, 
each existing particle (i) dies with probability a At, (ii) 
branches into two offspring with probability b At and (iii) 
diffuses, with diffusion constant D, with the remaining 
probability I — (a -I- 6) At. A schematic trajectory of the 
process is shown in Fig. ^ where Amax and Amin denote 
respectively the maximal displacements of the process up 
to time t in the positive and the negative direction. It is 
convenient to define these observables in their dimension¬ 
less forms aimax = Amax/\/-D/& and aimin = Amin/ \/D/b. 
Since the particle starts at the origin, a^max > 0 necessar- 
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ily and similarly ccmin < 0 necessarily. In the subsequent 
discussions we find it convenient to consider the positive 
quantities ccmax and — a^min as our basic random variables. 

As mentioned in the introduction, the marginal PDF 
of Xmax (and consequently that of —Xmin), has been ex¬ 
tensively studied for all a and b ii0- While in the 
supercritical phase (b > ^^this marginal PDF remains 
time dependent for all t [a, l7| , for 6 < a, it approaches a 
stationary form Pmarg{a;) which is known explicitly. It is 
convenient to express it in terms of its cumulative distri¬ 
bution 7Z(x) = PniaTAx')dx'. We set A = a/b— I. In 
the critical case A = 0 [ 3 ], 


n{x) = 




2 ■ 


( 1 ) 


Consequently, Pmaig{x) = —dTZ(x)/dx has the asymp¬ 
totic behaviors 


Pmarg(^) 


Pmarg(O) = \/3 : a; 0 


( 2 ) 


12 


, cc —>■ 00 . 


In the subcritical case A > 0 Q: 

^ • V 

-a; -I- sinh 


'R-{x) = -^csch' 


3A , 2 ( ,-1 



( 3 ) 


This result can further be simplified to give 


7^(x) 


a —1 

[—I -I- a cosh('\/Aa:) + \/— 1) sinh (-x/Aa;)] 


( 4 ) 


where a = 1 + 3A. 

Consequently, Pmarg(a:) = —dTZ{x)/dx has the asymp¬ 
totic behaviors 


{ Pmarg(O) — y -l-A, a;—>-0 

( 5 ) 

6Afe-2>*“h“iv^exp(-yAa;) , x ^ oo . 


By symmetry, — Xmin has the same marginal PDF 
Pmarg(a:) for A > 0. While Pmarg(a:) is thus well known, in 
this paper we compute the joint stationary PDF p{x, y) 
of X = Xmax and y = —Xmin for A > 0. One of our main 
results is to highlight the nonzero correlation between 
a^max and —Xmin cven in the stationary state. Indeed we 
show that 

P{x, y) ^ Pu ncorr [x,y) =Pmarg(a:)p marg (y) • (6) 

The stationary PDF of the dimensionless span C = 
s/^b = a;max — a;min can be computed from the joint 


PDF p(x, y) via the relation 

pOO pCO 

p(C) = / dx dyp{x,y)6{x + y-C) ■ (7) 

ao Jo 

We compute p{Q exactly for all A > 0 and find the 
following asymptotic behaviors. 

Critical case (A = 0): In this case we find 


piO 


A r- 

al = 87ry3 = 43.53118... 

C'^ 


( 8 ) 


Subcritical case (A > 0): Here we get 

P(C) ^ \ (9) 

exp(-\/AC) , C-too. 


where A = 12 A 


^3^-b \/l-f 3A/2 


-2 


Signatures of the correlation between Xmax and 
Xmin- Interestingly, one can show that these asymptotic 
behaviors of p{() for the critical ([8]) and the subcritical 
cases ([9]) carry the signatures of the correlation between 
a^max and Xmin (see also Figs. |8] and |9] below). In order to 
demonstrate this, we compute the asymptotic behaviors 
of p ((C) in the hypothetical case where one assumes that 
Xmax and Xmin are completely uncorrelated. Given that 
Xmax and —Xmin have the same PDF Pmargix) [obtained 
from Eq. ([IJ for A = 0 and from Eq. Q for A > 0], 
the span PDF Puncorr(C)! assuming that Xmax and Xmin 
are uncorrelated can be obtained by inserting p(x, y) = 
Puncorr{x,y) iuto Eq. dll) and is given by 


Puncorr(C) — f Pmarg ((^)Pmarg (C x'jdx . 
Jo 


For small C,, it behaves as 

Puncorr 


(C) ~ yLrg(O) C , when C ^ 0 . 


( 10 ) 


( 11 ) 


Substituting Pmarg(O) = y^2/3 + A from Eq. dS]) in Eq. 
dm gives 



C , when C —)■ 0 . 


( 12 ) 


Comparing this result with the exact one in Eqs. m and 
we see that, while both of them grow linearly for 
small (C, the slopes are different, reflecting the fact that 
aimax and Xmin are actually correlated. 

To investigate the large C behavior of Puncorr(C) in Eh- 
dm, we need to treat separately the critical (A = 0) 
and the subcritical (A > 0) cases - see Eqs. dll) and dl]). 
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In the critical case (A = 0), substituting the asymptotic 
behavior Pi„arg(a^) 12/x^ from Eq. (I5|) in Eq. (ITUl) . one 
gets for large ( 

24 

Puncorr(C) ^ TT ^ for A = 0 . (13) 

While this uncorrelated assumption correctly reproduces 
the decay ([8]), the prefactor 24 is different from the 
exact value A = Stt'/E = 43.53118... in Eq. again 
reflecting the nonzero correlation between Xmax and ccmin- 
On the other hand, for A > 0, one obtains from Eqs. ® 
and (fTU)) : 

Puncorr(C) ~ ^-4^ C exp (^-VA , for A > 0 . (14) 

Here also, the assumption of vanishing correlation cor¬ 
rectly reproduces the C"dependence oc i^exp c) 

of the right tail ([9]) but the amplitude is incorrect by a 
factor A/2, reflecting once again the presence of finite 
correlations between Xmax and ccmin- 


III. JOINT DISTRIBUTION OF THE 
MAXIMUM AND MINIMUM 

We are interested in the spatial extent of the BBM 
process up to time t. The process begins with a single 
particle at a; = 0 at time t = 0. We recall that the span 
of the process up to t, characterizing the spatial extent, 
is defined as s = Amax —-^min 7 where A„iax and Amin are 
respectively the maximum and minimum displacements 
of the process up to time t (see Fig. [T]). 

We start by defining the joint cumulative probability 

Q(A, Y, t) = Pr {Amax < A, Amin > -F; up to time t} . 

This has the simple interpretation as the probability that 
the process is confined within the box [—F, A] up to time 
t (see Fig. [21). The marginal cumulative distribution of 
the maximum can be obtained by taking the F —>■ oo 
limit of Q{X,Y,t)- Similarly the marginal cumulative 
distribution of the minimum is obtained by taking the 
A —>■ oo limit. The joint PDF P(A, F, t) of Amax = A 
and —Amin = F is then given by 

P{X,Y,t) = ^^QiX,Y,t) (15) 

where, by definition, A > 0 and F > 0. The PDF of the 
span s = A -f F is then given by 

pOO poo 

P{s,t)= / / dXdYS{X + Y - s)P{X,Y,t). (16) 

Jo Jo 


A. Backward Fokker-Planck equation for Q{X,Y,t) 

We derive a backward Fokker-Planck (BFP) equa¬ 
tion for Q(X,Y,t), following similar steps as in 
Refs. [3, O. We investigate how Q{X,Y,t) evolves 
into <5(A, Y,t + At). The goal is to derive a differential 
equation for the evolution of Q{X, Y, t). For this purpose 
we split the time interval [0, t + At] into two subintervals: 
[0, At] and [At,t + At]. We then take into account all 
possible stochastic events that take place in the first 
subinterval [0, At]. In [0, At], the particle at a; = 0 can 

A) split into two particles with probability bAt, result¬ 
ing in two BBM processes that are both confined within 
[—F, A] up to time t -|- At with probability Q^(A, F, t). 
The contribution from this term to Q{X,Y,t + At) is 
then bAt Q^(A, F, t). 

B) die with a probability aAt, leading to no particles 
at subsequent times. This event automatically ensures 
with probability 1 that the process remains confined 
within [—F, A] up to t -I- At. Hence it contributes a term 
aAt X 1 to Q(A, Y,t + At). 

C) diffuse with probability 1 — (5 -|- a)At, moving 
a distance Ax = T]{0)At in the first time step. This 
shifts the process by a distance Ax at the first time 
step. The probability that the resulting process is 
confined within [—F, A] up to time t -f At is then given 
by {Q{X — Ax, F -I- Ax,t))^(o). By the subscript r]{0) 
we denote an averaging over all possible values of the 
diffusive jump at the first time step. Hence this term 
contributes (1 — (5 -I- a) At) x {Q{X — Ax, Y + Ax, t)),,(o) 
to the final probability Q{X, Y,t + At). 

Adding the contributions from these three terms A), 
B) and C), we have 

Q{X,Y,t + At) = bAt Q^{X,Y,t) + aAt 
-f (1 — (6 -f a) At) {Q{X — Ax, Y + Ax, t))^(o))(17) 

where ry(t) is a Gaussian white noise process with the 
properties 

ivit)) = 0 , {vit)v<j')) = 2F'(5(t - t'). ( 18 ) 

Taylor expanding Eq. HU), using the properties of the 
noise in Eq. dl), and taking the limit At —> 0, we arrive 
at the exact BFP evolution equation 

-{b + a) Q{X,Y,t) + bQ^{X,Y,t). (19) 

Since at time t = 0, both the maximum and minimum of 
the process is at x = 0, the initial condition is 

Q(A,F,O) = 0(A)0(F) , (20) 
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where 0 is the Heaviside step function defined as 
1 for X > 0, 


©(x) = 


0 


for 


a; < 0 . 


( 21 ) 


At any time t > 0, the maximum ^rnax ^ 0 and the 
minimum —Amin > 0, leading to the boundary conditions 


Q(X,Y,t) = 


for 

for 


A < 0, 

y < 0 . 


It is actually convenient to work with 

R{X,Y,t) = l-Q{X,Y,t) 


( 22 ) 


(23) 


which denotes the complementary probability that the 
maximum or minimum up to time t is not within [—A, A]. 
Inserting Eq. (1^ into Eq. m we have 


(24) 

(25) 


dR{X,Y,t) fd dV 
dt ^\dX dYj 

+ {b-a) R{X,Y,t) - b R^{X,Y,t), 

with the initial conditions 

i?(A,r,o) = 0(-A)0(-y), 

and the boundary conditions 

' I for A < 0, 


R{X,Y,t) = 


I 


for 


y < 0. 


(26) 


B. Dimensionless Variables 


It is natural to consider the evolution equations in 
terms of dimensionless variables as follows 


A 


X = 


y = 


VWb’ 

Y 

~7Wb' 

T = 




(27) 


Similarly, we can define the dimensionless span of the 
process as 


C = x + y = 


.JWb 


(28) 


Our goal in this paper is to derive the stationary joint 
PDF of X and y and also the stationary PDF of C- In or¬ 
der to avoid a proliferation of symbols, we keep the same 
notation for the PDF’s of the unsealed and scaled vari¬ 
ables, with P(A, y, t) —>■ P(x, y, r) and P{s, t) —>• P(^, r). 


Similarly we have -R(A, A, t) —>■ i?(x, y, r). The distribu¬ 
tions of the scaled variables are related to the unsealed 
distributions as 


P(C, r) = P ( C = , T = bt] , 

/d76 /d 76’ y 

\ ( X Y 

P{x, y, t) = -T^P X = , y = . — T = bt\ . 




(29) 


In terms of these scaled variables Eq. (l24l) takes the 
simpler form 


dR{x,y,T) 

dr 


V clx dy ) 
-AR{x,y,T 


R{x,y,T) 

) - R^{x,y,T). 


(30) 


Eq. (1501) is a non-linear equation whose explicit solution 
at finite time t is hard to obtain analytically. In the 
supercritical case b > a, we expect this solution to be 
time dependent at all times t. However, for b < a, we 
show below that as r —>• oo, Eq. dini) admits a stationary 
solution for R{x, y, r) that can be computed explicitly. 
Using this solution, and Eqs. (HSl), (1211) and (1^ . the 
joint PDF of X and y can then be expressed as 

P{x, y, r) = y, r). (31) 

Finally, this joint PDF can be used to evaluate the PDF 
of the dimensionless span of this process defined in Eq. 
(HU), which is then given by 


P{C,r) 



dxdy5{x + y- C)P{x, y, r). 


(32) 


At large times, this PDF converges to the stationary dis¬ 
tribution p(C) = P{C,,T —>■ oo), which we analyse in detail 
in section 0 


IV. STATIONARY JOINT DISTRIBUTION OF 
Xmax AN^D Xmin FOR b ^ CL 

We now focus on the case b < a where the joint distri¬ 
bution R{x, y, r) in Eq. dliH) is expected to approach a 
stationary limit as r —>■ oo: 

TZ{x, y) = R{x, y,T ^ oo). (33) 

Setting the left hand side (Ihs) of Eq. (l30l) to 0 in the 
stationary limit gives 

( d d 

for X > 0 ,y > 0. Next, it is convenient to make a 
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FIG. 3: The change of variables {x, y} —>■ w}. 



change of variables 


( = x + y, 

v = x-y, (35) 

with C, G [0,oo) and v G [—CiC] (see Fig. [3]). Note that 
the variable C represents the dimensionless span of the 
process. In terms of these new variables Eq. (|M1) be¬ 
comes 

4(^) nic,v)-Anic,v)-n^Cv) = o, (36) 

valid in the regime v G [—C, +C] and C G [0, +oo) (see 
Fig. [3]). When v —>■ +C, i.e., y —0, the boundary con¬ 
dition R(x,y = 0,t) = 1 given in Eq. (1331) . translates 
into TZ((,() — 1. Similarly, when v —>■ —C, i.e., a: —>■ 0, 
the boundary condition R(x = 0,y,t) = 1 translates into 
7l(C, —() = 1. In addition, the solution must be symmet¬ 
ric around ?; = 0 (corresponding to x = y). Since v) 
is a cumulative probability, 0 < TZ{(,v) < 1. Conse¬ 
quently, for a fixed C, as v decreases from we expect 
that TZ{(,v) should decrease from its value TZ{(^,() = 1- 
By symmetry, as v increases from TZ{C,,v) should 
decrease from its value TZ{(,—Q = 1. Thus we ex¬ 
pect TZ{(,v), as a function of v for fixed C, is a smooth 
non-monotonic function, symmetric around u = 0 in 
~C ^ V < -bC, and with a minimum at n = 0 (see Fig. 
a. Assuming analyticity around the minimum at u = 0 
gives the condition 


a7^(c,u) 


dv 


= 0 . 




(37) 


Once we hnd the solution TZ{(^,v) of Eq. (1331) . using 
Eqs. (1311) and (1351) the joint PDF of ( and v can be 
expressed as 


19 9^, , 



^ncv). 

(38) 


FIG. 4: TZ{C, v) as a function of u G [—C, +C1 for different 
values of A For hxed C a) is a smooth non-monotonic 
function, symmetric around u = 0 in —^ < v < -bC and has 
a minimum at w = 0. The data plotted above corresponds to 
the case A = 0 and was obtained by numerically evaluating 
Eqs. (l5lTl and (l52l) . 


where the factor 1/2 in (I38p comes from the Jacobian of 
the transformation {x, y} —>■ {C, "c} (1351) such that 



f+C 

J ^ dvp{C„v) = I . 


(39) 


Fortunately, Eq. (|33)) can be integrated with respect 
to V upon multiplying by a factor , yielding 




^7^2(C,u) + i7^3(C,^;)+K;(C), (40) 


where k(/') is a yet unknown integration constant. To fix 
«;(/■), we use the condition in Eq. (1371) and arrive at 




^ {n\Cv)-n^{C0)) 

+ i (7^3(C,^;)-7^3(C,0)) . (41) 


Since the solution TZ((^,v), is symmetric about the 
?; = 0 line, it is sufficient to solve Eq. (HIT) for only the 
region v G [0,-bC] (or alternatively for v G [—C, 0]). We 
restrict on n G [0,-bC] where dTZ{(,v)/dv > 0 (see Fig. 
H. Taking the square root of Eq. HID and integrating, 
we obtain for u G [0, -bC]: 

_dr_ 

(^2_7e2(C,0)) + i(r3-7e3(C,0)) 

This equation can be conveniently expressed as 




92 92 

dv"^ 


( 42 ) 
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FIG. 5: The function ( 7 ( 7 , 2 ) for different values of 7 . For 
large 2 , Q{'^,z) saturates to a 7 dependent constant C( 7 ). 
The case 7 = 0 corresponds to the function C7(0, 2 ) analyzed 
in the critical case. The limiting behaviors are (7(0, 2 ) —>■ 0 as 

2—^1 and 0(0, 2 ) —>• C* = -^ « 2.4286 as 2 ^ 00 . 


y^(^^ V7^(C,0)’7^(C,0)y' 76’ ^ ^ 

where the bivariate function Q is defined by the integral 

z)= f /. 3 

Jl -1)+^ (x^ - 1) 

The above function 0 ( 7 , 2 ) can then be expressed as (us¬ 
ing the identity 3.138 of Ref. 0 ): 


0 ( 7 , 2 ) = 

_ \ _F 

(3 -b 27 ) 1/4 


2 tan 


2-1 1 


3-1-7 

73-b27 ’2V" 3-b27 


2 - 


where 2 > 1 , 7 > 0 and 

F(</>,fc)= [ 
Jo 


de 


0 \/l — k'^ sin^ 9 


(45) 


(46) 


is the elliptic integral of the first kind. In Fig. O we plot 
the function 0 ( 7 , 2 ) as a function of 2 for different values 
of 7 . Next, inserting the boundary condition 72.(7 =*=0 = 
1 in Eq. (PI) we have 


1 ^ ( 3A/2 1 

7^(7^^U(C,0)’72(C,0) 

This is an implicit equation for 72(C,0), the solution of 
which can then be injected in Eq. (P)) to solve for 72(71’) 



for all C and v. 


Critical Point. The computations become slightly more 
explicit exactly at the critical point a = b, i.e., A = 0. 
In this case, putting 7 = (3A/2)/72(7 0 ) = 0 in Eq. (ITTl) 
gives 


1 


=0 0 , 


72(77 


772(C,0)"' V''’^(C,o) 
where, from Eq. (HU) . 0 ( 0 , 2 ) is given by 
dx 


V 

7i’ 


0 ( 0,7 = 


j: 


'Jx'^ — 1 


2v0Fr(|) 2 


r(|) 


(48) 


117 1 




6 ’ 2 ’ 6 ’ 23 


where 2 F 1 is the usual Gauss hypergeometric function 
and 2 > 1. This function has the following asymptotic 
behaviors 


0 ( 0,7 


7^r(|) 2 

, 3 r(|) 


for 2 ^ 1 , 


for 2 —>■ 00 . 


(50) 


Similarly, putting A = 0 in Eq. (H71) determines 72(7 0) 
implicitly as the solution of 


=0 0 , 


772(C,0)"' V''’^(C,o) 
Dividing Eq. dH]) by m gives 

nc,v)\ v^f.. 1 


76' 


0 0 , 


72(70) 


V’ 72(70) 


(51) 


(52) 


where 0 ( 0 , 2 ) is given explicitly in Eq. (HHl) . The so¬ 
lution of Eq. (|52[) thus determines the cumulative joint 
distribution 72(7'(') in the critical regime. 


V. STATIONARY DISTRIBUTION OF THE 
SPAN 

The solution 72(7 '^) obtained from Eq. (H51) can next 
be used to compute the stationary joint PDF p(7 v) from 
Eq. (j38p . The stationary PDF of the dimensionless span 
7 denoted by p(C), is then obtained, for all A > 0 , by 
integrating over v as: 


r+C 

P{C) = J p{(,v)dv . (53) 

It is then easy to extract p(C), for all A > 0 , numeri¬ 
cally exactly from Eqs. (H51) . (H^ . (P)) and (1551) . As an 
example, we plot this PDF p{Q as a function of C for the 






















































FIG. 6: The function TZ{C,, 0) versus ^ in the critical regime 
derived using Eq. (EU, showing the limiting behaviors 
0 ) —>■ 1 as C —>• 0 and 7?.(C, 0 ) —^ as ^ ^ oo (dashed 

line) as predicted in Eq. (I60II . B « 35.3901 is defined in Eq. 
(I61II . Inset: Plot of 1 — 'n{C,, 0) showing the limiting behavior 
1 — ??.(((, 0 ) ~ (I) as ^ ^ 0 (dashed line) as predicted in 

Eq. ([5^. 


FIG. 7: Asymptotic behavior of TZ{C,,Q) in the subcritical 
regime. The plot shows 7?.(^, 0) for different values of A 
derived using Eq. (|i71). The dashed lines representing the 
asymptotic exponential behavior 7?.(^, 0)~Aexp as 

^ oo derived in Eq. (ED are indistinguishable from the 
theoretically obtained curves as they match exactly. Inset: 
Plot of 1 — 7?.(((, 0) showing the limiting behavior 1 — 0) ~ 

I (1 + A) as ^ 0 (dashed lines) as predicted in Eq. (1^ . 


critical case (A = 0) in Fig. |8]and for the subcritical case 
(A = 1) in Fig.O From the numerical plots for differ¬ 
ent values of A, one finds that for small C, p{C) increases 
linearly with a slope that depends on A. In contrast, for 
large p{Q has an algebraic tail p(C) oc for A = 0 
while it has an exponential tail for A > 0. In the next 
two subsections, we show that these asymptotic behav¬ 
iors of p(C), both for small and large (, can actually be 
extracted analytically for all A > 0. 


A. Asymptotic behavior of p{Q for C( 0 


in Eq. (l43ll with both C and v small gives 

nc, v) = l-^{l + A){e- v^) + 0(C^ V*). (57) 

Next, using Eq. (1551) . the joint PDE p{C,v) is given by 

= I(l+A)+OK’,i.*). (58) 

Substituting this expression in Eq. (l53ll gives 


As C —t 0 we have 7?.(C, 0) —>■ 1. Therefore we write 


7^(C,0) = l-e(C), (54) 


where e(C) is small. Substituting this in Eq. (lT7l) and 
expanding Q('y,z) in Eq. (HU) around z = I, we obtain 
(using the notation e = e{Q) 


2\/e 

y/3(l -f A) 




ye' 


Therefore to leading order in (J) we have 


(55) 


p(C) = i(l + A)C + 0(y) , (59) 

which yields the small C behavior announced in Eq. m, 
for A = 0, and in Eq. (E|) for A > 0. The asymptotic 
linear growth for small C is shown, for the critical case 
(A = 0), in the inset of Eig. [8] and for the subcritical 
case (for A = 0.1,1,10) in the inset of Eig. Oa.. As 
discussed in section m the amplitude of this linear term 
in (1591) carries the signature of the correlation between 

^max and Xniin- 


7^(c,o) = i-i(i+A)y+o(y). (56) 


B. Asymptotic behavior of p(c() for ^ oo 


This limiting behavior for the critical case (A = 0) is 
illustrated in the inset of Eig. [SI and for the subcritical 
case in the inset of Fig. [T] Performing the same analysis 


In this subsection we extract analytically the large 
tails of p{() both for the critical (A = 0) as well as for 
the subcritical case (A > 0). 
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1. Critical point (A — 0) 

We start by analyzing Eq. eiD in the limit ( —i oo. 
In this limit 7l((, 0) is small. We then need to analyze 
Q(0,z) for large z. Inserting the asymptotic behavior of 
Q(0,z) in Eq. (1501) into Eq (I5T]) . we obtain 

7^(C,0) = ^+O(^l), (60) 

where 

S = 6 Ri 35.3901, with C* = — (61) 

3 r(§) ^ ^ 

This asymptotic behavior of 7l((,0) is illustrated in 
Fig. |6l Having thus determined 7l(C, 0) for large (, we 
now investigate 7Z((, v) for large (. Our aim is to extract 
p(C) for large ( from Eqs. (|05)) and (|50)) . We note that 
Eq. (1531) involves an integral over v and this integral 
is dominated by u ^ c). Hence we need to investigate 
TZ{C,v) in the scaling limit C oo, u —?> oo but keeping 
C/t) fixed. 

When (( —>• oo, TZ{C, 0) —>■ 0 as in Eq. (1001) . As a result, 
the argument l/7^(^, 0 ) of 1 /( 0 , l/7?.(^, 0 )) on the right 
hand side of Eq. (|52l) goes to oo. From Eq. (|501) . we see 
that 1/(0,z —>■ oo) = C* where C* is given in Eq. (10T1) . 
Hence, in the scaling limit C —>■ oo, n —>■ oo keeping (/v 
fixed, Eq. (l52l) becomes 



FIG. 8: Theoretical stationary PDF of the dimensionless span 
p{() (solid line) in the critical regime (A = 0) calculated by 
numerically evaluating Eqs. dSIl), dSH), IMI) and ((531). The 
(blue) dotted line corresponds to a plot of puncorr(C) evaluated 
from Eqs. m and uni which is the result for the span dis¬ 
tribution obtained under the assumption that Xmax and Xmin 
are uncorrelated. The Inset displays the same data in log-log 
scale showing the two limiting behaviors (dashed lines), linear 
P{0 ~ (1) C for small spans, and p{Q ~ for large spans 
(El), with A = SnVS ~ 43.53118 .... 


gives finally 


g 



nc,v) \ 

7^(C,0)y' 



Inverting the above Eq. we get 




cC* 


dz 


6J^iz)+QzJ^' (z) + 

(z^ - _ (-66) 


nc,v) 

7^(C,0) 



(63) 


where J^{z) is defined as the inverse function of C/( 0 ,z). 
Substituting TZ{(, 0) ^ B/('^ from Eq. (1501) gives the final 
scaling limit expression of the joint cumulative distribu¬ 
tion 


TZ{Cv) = ^J^(c*^y (64) 

Inserting this expression into Eq. (1381) we arrive at the 
joint PDF 


Hence, we obtain 

A 

p{0^ for large C, (67) 

with 


rC” 


A = -6C* 


dz 


6A{z) + 6zA'{z) -b (z^ - C*‘^)A”{z) 


( 68 ) 

It turns out that the integral in Eq. (1051) can be per¬ 
formed explicitly (see Appendix!^ and the final answer 
for the amplitude is amazingly simple 


^ = 871^3 = 43.53118.... (69) 


B 


p(C,v) = -- 


2 LC 


-tA C*- +6C*—A' C*- 


+C* 




■ (65) 


Substituting this expression for p(C, v) in Eq. (1551) 


Thus the leading asymptotic behavior for large () is 


p{C) 


877^3 


(70) 


as announced in Eq. This large () behavior of the 
span PDF is shown in Fig. |51 where we plot the numer- 
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ically exact distribution p{() (extracted from Eqs. (ISTI) . 
(EH), El and (ESI), along with the asymptotic power 
law tail derived analytically in Eq. EH- At large 
the exact distribution shows a good agreement with the 
asymptotic behaviour. We also display Puncorr(C): the 
span PDF obtained assuming that ccmax and Xmm are un¬ 
correlated, which decays with the same power, but with 
a different prefactor provided in Eq. m- The ampli¬ 
tude A = S'K'/Z is thus nontrivial due to the remnant 
non-vanishing correlation between ccmax and a^min in the 
stationary state (see also the discussion at the end of 
section El. 


2. Subcritical case (A> Q) 

In the subcritical regime, the large C analysis can be 
done along the same lines as before for the critical case 
but the analysis is a bit more complicated and most of 
the details have been relegated to Appendix [BJ 

As in the critical case, here also we first need to com¬ 
pute the large C, behavior of 72.(C, 0) from Eq. (1471) with 
A > 0. This is done in Appendix [BJ We find (see Eq. 

dEH) 


7e(C,0) - Aexp (-^C I , as C oo 


where A = 


12A 


vW2 + \/1 + 3A/2 


2 • 


(71) 


In Fig. Ewe show a plot of 7^(C,0) obtained by numer¬ 
ically evaluating Eq. (iTTl) for different values of A. For 
large C, this shows a good agreement with the asymp¬ 
totic behavior in Eq. ED- Using this asymptotic tail 
of 7^(C, 0) from Eq. (ITT]) , we then analyse the asymp¬ 
totic behavior of TZ{(,v) in Eq. (H51) . in the scaling limit 
C —>■ oo, u —>■ oo, keeping (/v fixed. We find 

7^(C, v) = HoiC, v) + 7^2(C, 0)i7i(C, v)+0 {TZ^iC, 0)) . 

(72) 

where both Ho{C,v) and Hi{(^,v) are functions of m = 
C — V only, whose expressions are provided in Eq. (IB16I) 
and Eq. (IB17I) respectively. From TZ{(,v), we can then 
obtain the joint PDF p{(, v) from Eq. (l38l) and eventu¬ 
ally p(C) from Eq. (1551) . Following this procedure, we 
finally obtain the large ( behavior of p{Q for A > 0 (see 
Appendix iBl for details): 

p{C) = ^Ce-'^Hl + OiC^)) , (73) 

with A given in Eq. (ITT]) , as announced in Eq. ED- Here 
again, as discussed in section El the amplitude A^/2 in 
Eq. (1731) bears the signatures of the correlations between 
a^max and In Fig. [5] we show a plot of p(C), for 

the subcritical case (for A = 1), obtained by numerically 
evaluating Eqs. (155|) - (H71) . (|5^ and (|55)) - For compari- 




FIG. 9: Theoretical stationary PDF of the dimensionless span 
p{() in the subcritical regime. In a. we display p{Q (solid 
line), for A = 1, calculated by numerically evaluating Eqs. 
(|43ll. (HzD, EHD and EH- The (blue) dotted line corresponds 
to puncorr(C) evaluated from Eqs. ED and (1101) which is the 
result for the span distribution obtained under the assump¬ 
tion that Xmax and Xmin are uncorrelated. The Inset displays 
the same data in log-log scale showing the approach to the 
0 linear behavior (dashed line) derived in Eq. (15911 . In 
b. we display the stationary span PDF p{Q and puncorr(C) for 
A = 1 in log-linear scale. The dashed line showing the ^ ^ oo 
asymptotic behavior ~ ^(^exp (—Q derived in Eq. dZS) 
is indistinguishable from the theoretically obtained curve as 
they match exactly. The Inset shows the same data high¬ 
lighting the difference between the asymptotic decays of p(^) 
and Puncorr(C) (the asymptotic decay of Puncorr(C) is provided 
in Eq. (fTil)'). 


son, we also show a plot of Puncorr(C) obtained from Eqs. 
ED and EHD which corresponds to the PDF of the span 
obtained by assuming that Xmax and XnPn are indepen¬ 
dent. We also display the agreement of the asymptotic 
behaviors for C 0 and (->-00 derived in Eqs. dMD and 
(1751) with these numerically exact PDFs. 
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VI. MONTE CARLO SIMULATIONS 

We have performed numerical simulations of the BBM 
and numerically computed the PDF of the span at dif¬ 
ferent times t. Directly simulating the BBM model is in 
general hard to do in the supercritical regime {b > a) 
where there is an exponential proliferation of particles. 
In this case one has to resort to numerically evaluating 
the non-linear FKPP-type equations to extract the be¬ 
havior of the PDF at large times However, in the 

critical (b = a) and subcritical (b < a) cases, it is pos¬ 
sible to obtain very good statistics by performing direct 
Monte Carlo simulations of the process la¬ 
in Fig. [TU] a. we show our numerical results at crit¬ 
icality (A = 0) for the PDF of the dimensionless span 
P((, T = bt) at large time t = 100 (note that the discrete 
time step was set to At = 0.0001 in our simulations) 
with D — 1 and a = b = I. These data show a very good 
agreement with our theoretical predictions for the sta¬ 
tionary PDF p{() extracted numerically from Eqs. ([Hi]) . 
(|5^ . (p8|) and ([^ . In Fig-lTOlb. we show the same quan¬ 
tity in the subcritical regime for A = I. Here again, we 
observe a very good agreement between the Monte Carlo 
simulations and our exact theoretical results extracted 
numerically from Eqs. (H51) . (HTl) . (1551) and (1551) . except 
in the small C region where discretization effects become 
important. We have checked that as we decrease the size 
of the time steps, the results from our simulations match 
more closely with our theoretical predictions. 

Finally, we have also numerically studied the finite 
time behavior of the solution of Eq. (1241) . In particu¬ 
lar, for A = 0, our data for finite times indicate that 
P{(,t) takes the scaling form P{(,t) ^ p{()S{C/^/ t), 
where S{u) is a rapidly decaying function as it —^ oo. The 
curvature of the PDF at large C in the inset of Fig. llUla. 
is a feature that emerges due to the finite time nature of 
our measurement, and our numerical solutions reproduce 
it exactly. It would be interesting to analytically study 
this finite time behavior which certainly deserves further 
investigation. 


VII. CONCLUSION 

In summary, we have obtained exact results for the 
stationary joint PDF p(a;, y) of the dimensionless max¬ 
imal displacements Xmax = ^max/ \/D/b and ccmin = 
X-caSxij\fl5]b up to time t for the one-dimensional BBM 
in the critical (a = b) and subcritical cases (a > b) (see 
Figs. [T]and[5|). In both cases we found that the corre¬ 
lation between ccmax and Xmin remain nonzero, even in 
the stationary state. From this joint PDF we have com¬ 
puted exactly the PDF p{() of the (dimensionless) span, 
C = — Xjnin)/\/D/b, which provides an estimate of 

the spatial extent of the process. We demonstrated that 
p(C) carries the signatures of the correlation between the 
two extreme displacements ccmax and Xmin, which can be 



0 5 10 15 20 

b. C 

FIG. 10: a. Probability distribution function of the dimen¬ 
sionless span p(^) extracted from Monte Carlo simulations 
(open circles) in the critical case (A = 0). Here t = 100, 
D = 1, a = fe = 1, and dt = 0.0001. The data is averaged over 
5 X 10^ realizations. The bold line represents the stationary 
theoretical PDF p{() extracted numerically from Eqs. (1511) . 
dSl), dMJ and dSSJ. The Inset displays the same data over 
a wider range in log-log scale. The curvatnre of the PDF for 
large arises due to the finite time nature of our measure¬ 
ment. 

b. Probability distribution function of the dimensionless span 
p{() extracted from Monte Carlo simulations (open circles) in 
the subcritical regime. Here t = 100, D — l,a = 2, b — 1 (i.e. 
A = 1), and dt = 0.0001. The data is averaged over 5 x 10^ 
realizations. The bold line represents the stationary theoreti¬ 
cal PDF p(C) extracted numerically from Eqs. (|43l) . (|47|l . (1381) 
and (1531) . The Inset displays the same data in log-linear scale 
showing the approach to the exponential behaviour for large 

as predicted in Eq. dZS]). 


seen for instance in the asymptotic behaviors ofp(C) both 
for small and large arguments diiii). 

The span is an interesting physical observable asso¬ 
ciated with BBM, which has several potential applica¬ 
tions, for example in the context of epidemic spreads [l3l| . 
Moreover, our results are also interesting from the general 
point of view of extreme value statistics (EVS) of strongly 
correlated variables. It was indeed recently demonstrated 
that random walks and Brownian motion (see e.g.. Refs. 
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for recent studies) are interesting laborato¬ 
ries to test the effects of correlations on EVS, beyond 
the well known case of independent and identical ran¬ 
dom variables [2^. In that respect, the results for the 
one-dimensional BBM obtained in the present paper con¬ 
stitute an interesting instance of a strongly correlated 
multi-particle system where the correlation between ex¬ 
treme values can be computed analytically. 

In this paper we have restricted ourselves to comput¬ 
ing the span distribution for the critical {b = a) and 
subcritical cases {b < a). The computation was feasible 
because the span distribution becomes stationary at late 
times t in these cases. In contrast, in the supercritical 
case (b > a) the span distribution will always be time 
dependent and it would be interesting to compute this 
distribution exactly. The recent developments @,i, mi¬ 
ls^ in the supercritical case may shed some light on this 
outstanding problem. 

It would also be interesting to extend these calculations 
to branching processes where the particles can split into 
k > 2 particles at each time step, which can be treated 
using the techniques developed in our paper. 


Acknowledgements 


J^{z) and its derivative can then be expressed as follows 
J-(z) = r, 

jz'(z) = -^4-7 = - 1- (A3) 

4?S{0,r) 

Similarly we can represent the integral in Eq. (IA2|) in 
terms of the following function. We define 


pC* poo 

/ dz' J^{z') = n{r)= / dr'- 
J z Jr 


+ (A4) 


^(i) 

r 


Inserting the above expressions into Eq. (IA2I) . we arrive 
at the following exact expression for the coefficient 


A = X{r) 


OO 

1 


(AS) 


where the integrand X{r) is defined as 


I(r) = 


- 6 C* 


2'H{r) + 4 rQ{Q, r) + (^1(0, r f - 


^3 _ 1 


(A6) 
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Appendix A: Critical Prefactor 


We next examine the limiting behaviors of the integrand 
X{r) in (IA6I) . Using the expressions of the functions 
G{0,r) in Eq. (l4^ and 'H{r) in Eq. (IA4I1 . we find 



(A7) 


In this appendix we derive an exact and simple expres¬ 
sion for the amplitude A of the power law decay of the 
span PDF in the critical regime, given in Eq. (1681) in the 
text. We first rewrite the integral in Eq. (1551) as 


A = 


- 6 C- 


L° ‘‘‘i 


2A{z)+ 4zT{z) + {z^ -C*‘^)A'{z) 


(Al) 


Integrating with respect to z we arrive at 


Using the above expressions we obtain the limiting be¬ 
haviors 

liml(r) —7> 0 and lim I{r) —>■ 87r-\/3 . 

r—>-1 r—>-oo 

Finally, inserting these into Eq. (|A5I) . we obtain the 
following exact value for the coefficient 

yl = 87ry3- 43.53118..., (A8) 

as announced in the text (l69l) . 


A=-6C* 


4zJ-(z) -f {z"^ -C*^)A'{z) 
c* 


-12 


C* f dzF{z). (A2) 
Jo 


Next, it is convenient to perform the change of variable 
r = T(z), with z = 1/(0,r) and r € [l,oo] when z € 
[0,C*], where the function G{0,r) is defined in Eq. (H51) . 


Appendix B: Asymptotic behaviors of the PDF of 
the span in the subcritical case A > 0 

We recall that the PDF of the span p{Q is given by 

Eq. dag): 

p{C) = j p{C,v)dv = 2j p{C,v)dv , (Bl) 
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where we have used p(C, v) = p{(, —v) and where p{(, v) 
is given by 


19 9 
2 dx dy' 


piC.v) = = ^ ( ^ - 7e(c,u). 

(B2) 

The function TZ{C,, v) is itself determined implicitly by 
Eq. (gSl) 


1 ^ ( 3A/2 'RJ^\ ^ 

U(C,o)’7e(c,o)y' ye’ 

where Gi'J, z) is given by 

z) = [ . 

Ji y(a:3 — 1) + 7 (a;^ — 1) 


(B3) 


(B4) 


The goal is to extract the behavior of the function TZ{C, v) 
from Eqs. (IB3I) . (IB4I1 in the limit of large ( and v with 
^ - as the integral over v in Eq. m is dominated 

by u -- C- 


First, it is convenient to rewrite the function C/( 7 , z) in 
(lB4)) as 


Gh,z) = Goo{i) - G{i,z) 


(B5) 


where 


^ 00 ( 7 ) = 
G{i,z) = 


dx 


1 y(a;3 - 1 ) + 7 (a ;2 - 1 ) 

dx 

y(a:3 - 1 ) + 7 (a ;2 - 1 ) 


(B 6 ) 


In the following, we will need the asymptotic behavior of 
1 / 00 ( 7 ) for large 7 : 


^;c^( 7 ) = ^(log7 + 31og2)-^ + (9(^) . 

(B7) 

On the other hand, we also need the asymptotic expan¬ 
sion of the function ^( 7 , z) in (IB 6 I) in the limit 'y ^ 00 , 
z —>■ 00 , keeping the ratio z = a fixed. This expan¬ 
sion can be obtained straightforwardly from the integral 
representation given in (IB 6 |) by performing the change of 
variable x = zu. One finds 


O 


2 sinh“^(Y^) 


where 


^i(a) = 


du 


1 ^/{y^ av?) \/cr 

In a -I- In 4 
\f0i 


(B 8 ) 


(B9) 


as a —>■ 00 


and 


G2{a) = 


du 


au 


20/2 


yS {2a^ — 5a — 15) -|- ISyo; -|- 1 sinh ^ {y/a) 

8ya®y7+T) 



as a —>■ 00 . 


(BIO) 


Solving Eq. (IB3I) for TZ{(, v) requires the knowledge of 
TZ{(, 0) which we first study, in the large C limit. This 
function 7?.(y 0) satisfies (as given in Eq. (07])) 

1 ^ ( 3A/2 1 

y^(yy^U(C,o)’7^(C,o) 

In the limit () —>• 00 , one expects that TZ{C,,0) —?> 0 
hence we need the asymptotic behavior of t/( 7 , z) for 
7 = (3A/2)/72.(C, 0) —>■ 00 , z = 1/7?.(C, 0) —>■ 00 keeping 
a = "fjz = 3A/2 fixed. Using Eq. (IB5I) together with 
the asymptotic expansion in Eq. (IB 8 I) at lowest order - 
i.e, retaining only Gi{a) - one finds from Eq. (IBIll) 


ye’ 


(Bii) 


7?.(y 0)~Aexp 



as C 00 


where A = 


12A 


vW2-byi + 3A/2 


(BI2) 


as given in Eq. dni in the text. 

We now study the asymptotic expansion of TZ{(,v), 
from Eq. (1531) and using the expansion of 7^(y0) ob¬ 
tained above (IB12I) . Here we need the asymptotic be¬ 
havior of C/( 7 , z) for 

7 = (3A/2)/7^(C, 0 ) ^ 00 , z = nc, v)/niC, 0 ) ^ 00 


keeping a = 7 / 2 : = 3A/{2TZ{(,v)) fixed. Inserting the 
asymptotic expansions obtained above (|B7ll - (IB8p in Eq. 
(IB3I) one finds 


I 


/ 3A \ nHc,0) ~ f 3A Y 

V27^(c,yy 7^yc,y V27e(c,yy 


2 y 2 

y^ 


sinh 



C-v 

ye 


(B13) 


which is valid up to terms of order 7?.^(0, y, which is small 
when (—>■ 00 [see Eq. (IBI2|) ]. Hence from Eq. (IB13I) . 
one expects that 7Z((, v) admits the following expansion, 
for C —>■ 00 , —>■ 00 with V ^ 

7^(C, v) = HoiC, v) + 7^yc, 0)i7i(C, v)+0 {TZ^C, 0 )) . 

(BI4) 

As we will see, to obtain the asymptotic behavior of p{() 
to lowest non-trivial order for large y we need to compute 
7?.(C,u) up to order 7?.^(y0), i.e. we need to compute 
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both Ho{(,v) and Hi{(,v) in Eq. (IB14I) . 

Computation of Ho{(,v). Neglecting the second term in 
the l.h.s. of Eq. (|B13I) . one obtains the following equation 
for Hq = Ho{C,v): 


1 




3A 


2^/2 


sinh-^ 


C-v 

V6 


(B15) 

which clearly shows that Hq is a function of C, — v, 
^o(Ci '^) = Ho{u = f — v) where Hq{u) is given by 


Ho{u) = 


3A 


2 sinh^ ^ 




(B16) 


which is obtained by injecting the explicit expression of 
Gi{a) in (IB9I1 into Eq. (IB15I1 . 


Computation of Hi(f,v). By inserting the expansion 
(IB14P into Eq. (IB13I) and expanding up to order 
one obtains that is also a function 

of u = f — V only, given by 


Hi{u) = :^Z{u) 


G2{Z{u)) 


with Z(u) = 


^^GiiZiu)) 

3A 


z{u)gi {z{u)) 


2Ho{u) 


(B17) 


where the functions t/i(a) and t/ 2 (<a) are given in Eq. 
(|B9I) and Eq. (IBlOjl respectively and Gi (a) = -£^Gi{a)- 


Computation of the PDF of the span p(f) for large f. 
From Eq. (IB2I) together with the expansion in Eq. (IB14I) , 
we obtain the joint PDF p{f, v) as 


PiC,v) = ^{dv - d^)n{C,v) = 

- (7^'(C, 0))" + 7^(C, 0)7^"(C,0)i7i(C - v) 

+27^(C,0)7^'(C,0)i^^(C-u))] +O(7^3(C,0)) , 

(B18) 


where we have used — d^)Ho{f,v) = 0 as Ho{(,v) = 
Hq(u = f — v) depends only on u = f — v. Using the 
large f expansion of TZ{(,0) in Eq. (jB12l) . one finds, 
from (IBISP 


p{C,v) 


(2VAH[{C -v)- Ai7i(C - u)) 
+o(e-^‘^^ . (B19) 


illustrated in Fig. 1111 where we plot this distribution ob- 



FIG. 11: Plot of the joint PDF p{Cv), obtained by numeri¬ 
cally evaluating Eqs. (|43ll. (|47ll. (|351l and (|53ll. as a function 
of C — w, normalized by for a fixed f (which we 

choose to be a typical large value ( = exp(3)). The dashed 
lines denoting the function 2y/AH[{(^ — v) — AHi{(^ — v) (de¬ 
rived in Eq. (IB19I) 1 for different A are indistinguishable from 
the theoretically obtained curves as they match exactly. This 
function tends to 1 as ^ — u —^ oo, as predicted by Eq. (IB21II . 


tained by numerically evaluating Eqs. dlSl), gzl) and 
(IMl) for a fixed C, (chosen to be a typical large value 
C = exp(3)). We find a very good agreement between 
the exact PDF and the asymptotic behavior derived in 
Eq. (IB19I) for all A. Finally, inserting this asymptotic 
behavior (jB19p of the joint PDF into Eq. (|Bll) yields the 
PDF of the span which is given by 

piO = 

.(B20) 

Using Eq. (IB16I) together with Eq. (IB17I) . one can show, 
for instance using Mathematica that 

lim (2VAH[{u) - AHi{u)^ = ^ . (B21) 

This limiting behaviour (IB21I) is illustrated in Fig. [TTlfor 
two different values A = 1 and A = 10. Finally, inserting 
this large C limiting behaviour (IB21I) into Eq. (IB20|) . we 
arrive at 


(2VAH[{u) - AHi{u)^ 


p{C) = ^Ce-^ ^ (1 + 0(r')) , (B22) 


This behaviour of the joint PDF p{f, v) for large ( is 


as given in Eq. (1751) in the text. 
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